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Abstract 

Complete one- loop results are obtained for the class of processes Xi ~^ ^j^^ 
MSSM where all parameters entering this process beyond lowest order are allowed to 
have arbitrary CP-violating phases. An on-shell renormalisation scheme is worked out 
for the chargino-neutralino sector that properly takes account of imaginary parts aris- 
ing from complex parameters and from absorptive parts of loop integrals. The genuine 
vertex contributions to the neutralino decay amplitudes are combined with two-loop 
propagator-type corrections for the outgoing Higgs boson. In this way the currently 
most precise prediction for this class of processes is obtained. The numerical impact of 
the genuine vertex corrections is studied for several examples of CP-conserving and CV- 
violating scenarios. We find that significant effects on the decay widths and branching 
ratios are possible even in the CP-conserving MSSM. In the CP-violating CPX bench- 
mark scenario the corrections to the decay width are found to be particularly large, 
namely, of order 45% for a Higgs mass of 40GeV. This parameter region of the CPX 
scenario where a very light Higgs boson is unexcluded by present data is analysed in 
detail. We find that in this parameter region, which will be difficult to cover by stan- 
dard Higgs search channels at the LHC, the branching ratio for the decay X2 ~^ Xi^i 
is large. This may offer good prospects to detect such a light Higgs boson in cascade 
decays of supersymmetric particles. 
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1 Introduction 



One of the main goals of physics at the Large Hadron CoUider (LHC) will be to shed light on 
the mechanism of electroweak symmetry breaking (EWSB). The most popular realisation 
of EWSB in theoretical models is the Higgs mechanism, which gives rise to at least one 
fundamental scalar particle in the spectrum. While in the Standard Model (SM), Higgs 
physics is determined by a single parameter, the mass of the Higgs boson, Mh, a much 
richer phenomenology is possible in extensions of the SM. 

In the Minimal Supcrsymmetric Standard Model (MSSM), the Higgs sector is charac- 
terised at lowest order by two new parameters instead of one. The spectrum contains five 
physical Higgs bosons, the properties of which may differ significantly from those of a SM 
Higgs. At lowest order the neutral Higgs bosons of the MSSM are CP-eigenstates, so that 
there are two CT'-even Higgs bosons, h and H, a CV-odd Higgs boson. A, and two charged 
Higgs bosons, H^. Higher-order contributions in the MSSM Higgs sector yield large cor- 
rections to the masses and couplings, and can also induce CP-violation, so that mixing can 
occur between h, H and A in the general case of complex SUSY-breaking parameters. If 
the mixing between the three neutral mass eigenstates, denoted hi, /i2 and h^, is such that 
the coupling of the hghtest Higgs boson, hi, to gauge bosons is significantly suppressed, this 
state can be very light without being in conflict with the exclusion bounds from the Higgs 
searches at LEP [1,2] and the Tevatron [3]. In particular, in the CPX benchmark scenario [4] 
an uncxcludcd region remains in which Af/i^ ~ 45 GcV and tan/? ~ 7 [2] (sec also Ref. [5] 
for a recent reevaluation with improved theoretical predictions). This unexcluded parameter 
region with a very light Higgs boson will also be difficult to cover at the LHC with the 
standard search channels [6-8]. 

While on the one hand a supcrsymmetric (SUSY) scenario such as the CPX scenario may 
have much worse prospects compared to the SM case for Higgs searches at the LHC in the 
standard channels, on the other hand additional Higgs production channels involving SUSY 
particles may occur in such a case. In cascade decays of heavier SUSY particles down to the 
lightest supcrsymmetric particle (LSP), Higgs bosons can in particular be produced in decays 
of neutralinos and charginos, via xf ~^ h, H or A and xf ~^ X^H^, see e.g. Refs. [9,10] 
for studies of these channels at the LHC in the MSSM with real parameters. These channels 
have also attracted recent interest for studies of scenarios with non-universal gaugino masses 
[11-13]. In the parameter regions of the CMSSM (the constrained MSSM) and the NUHM 
(a generalisation of the CMSSM with a non-universal Higgs mass parameter) which are 
currently favoured by electroweak precision data, i?-physics observables and cosmological 
data, an early discovery of the light Higgs boson from a neutralino decay in a SUSY cascade 
could be possible [14]. A related process to the production of a Higgs boson in the decay 
of a neutralino is the decay of a heavy Higgs boson into two neutrahnos, H,A ^ XiX°- 
This process, with a possible signature of four leptons plus missing energy, can also be 
phenomenologically important [15,16]. Concerning theoretical predictions for this class of 
processes, partial one-loop results have been published previously for the decays H,A^ "XaX^ 
in both the Feynman-diagrammatic [17, 18] and effective potential [19] approaches. These 
predictions did not include the full MSSM, and the Feynman-diagrammatic calculations were 
restricted to the case of real parameters. 

In the present paper we obtain predictions for decays of a heavier neutralino into a lighter 
neutralino and a neutral Higgs boson in the MSSM with complex parameters, i.e. we con- 
sider the class of processes x° —>■ x'jha, where ha = hi,h2,h3. Our calculations are also 
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applicable to the related class of processes ha XiX^- Since higher-order contributions in 
the MSSM Higgs sector are known to be large, a proper inclusion of Higgs-sector corrections 
is indispensable for a reliable prediction of this class of processes. The process-independent 
corrections to the mass of the outgoing Higgs boson and to the Higgs wave function normal- 
isation can be incorporated via an effective Born-type prediction for the neutralino decay 
process, see Refs. [5,20,21]. The genuine (process-specific) vertex corrections can also be 
very important. This has recently been demonstrated in Ref. [5] for Higgs cascade decay 
processes, ha —>■ hbhc, in the CPX scenario, where the genuine vertex corrections were found 
to give rise to drastic changes in the decay widths compared to the effective Born-type pre- 
dictions. In the neutralino decay processes, comprising just one instead of three external 
Higgs bosons, the genuine vertex corrections are not expected to be quite as large as for the 
Higgs cascade decays, but their effects can nevertheless be expected to be non-negligible. 

We use the Feynman-diagrammatic approach to evaluate higher-order contributions to 
the processes x^* X^ha- Specifically, we compute the vertex corrections at the one-loop 
level, taking into account the contributions from all MSSM particles, and we combine these 
results with state-of-the-art two-loop propagator-type corrections as implemented in the 
code FeynHiggs [20, 22-24] . In this way the currently most precise prediction for this class 
of processes is obtained. We focus our treatment of CP-violating phases on those that are 
most relevant for Higgs phenomenology, namely the phases of the trilinear couplings of the 
third generation, (f)At b t j ^^^^ ^^c gluino phase, (these are also the phases chosen to 
be non-zero in the CPX benchmark scenario [4]; the gluino phase enters the predictions 
for the neutralino decays via two- loop Higgs propagator-type contributions). We use an 
on-shell scheme for the renormalisation in the chargino-neutralino sector. In the MSSM 
with complex parameters care has to be taken in the treatment of absorptive parts of loop 
integrals and imaginary parts of MSSM parameters, since products of such contributions can 
enter predictions for physical observables already at the one-loop level. We have worked out 
a scheme for the renormalisation in the chargino-neutralino sector where in- and outgoing 
fermions receive different field renormalisation constants. 

In our numerical discussion we concentrate in particular on the parameter region in the 
CPX benchmark scenario where a light Higgs boson is unexcluded by current data (see also 
Refs. [25-27] for discussions of other possible LHC search channels to access this parameter 
region), but we also give examples for the CP-conserving case. Based on our results, we 
investigate the phenomenology of Higgs searches at the LHC in the channels x° x'jha- 
We briefiy discuss the prospects for covering the unexcluded parameter region of the CPX 
scenario in this way. 

2 Lowest-order Result, Notations and Conventions 

We first lay out our notation for the Higgs and chargino-neutralino sectors of the MSSM, 
and use this to write down a formula for the tree-level decay width for x° — > x^h^, where 
h^ is one of the neutral MSSM Higgs bosons, h, H or A. We also include notation for the 
sfermion sector which enters the process at the one-loop level. 
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2.1 Higgs Sector 

In the Higgs sector we follow the conventions of Refs. [5, 20]. We write the two Higgs doublets 
at tree level as 

V -01 y ' V ^2+^(02 + ^X2) J ^ ' 

The tree-level physical states h, H, A, and unphysical Goldstone states G, are ob- 
tained from rotations by the mixing angles a, /3„ and (5c as shown, 
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As indicated by the null entries in the 4x4 mixing matrix above, at tree level there is no CP- 
violating mixing between the neutral Higgs bosons. Minimization of the Higgs potential and 
the requirement of vanishing tadpoles at tree level renders the phase = and (3n = Pc = (3, 
where tan/5 = V2/V1 is the ratio of the Higgs vacuum expectation values. The Higgs sector 
is characterised by two input parameters (besides the gauge couplings) , conveniently chosen 
as tan /3 and one of the Higgs-boson masses. For the latter, the most convenient choice in 
the case where the SUSY-breaking parameters are allowed to be complex is the mass of the 
charged Higgs boson, Mh±, since the three neutral Higgs bosons mix with each other once 
higher-order corrections are taken into account. 



2.2 Chargino and Neutralino Sector 

At tree level, the physical chargino states, xf (i = 1,2), are Dirac spinors constructed from 
the mass eigenstates of the 2x2 complex mass matrix X, which reads, in the wino-higgsino 
basis, 

_ f M2 V2MwsmP 
~ \ V2MwC0sP II 

where M2 and fi are the wino and higgsino mass parameters, respectively. The off-diagonal 
elements depend on parameters from other sectors, namely tan/5 and Mw, the mass of the 
W boson. The mass matrix is diagonalised by two 2x2 complex unitary matrices U and V, 
where U*XV^ — diag(m-±, m-±). Similarly, the neutralinos (i = 1,2,3,4) are Majorana 
spinors constructed from mass eigenstates of the 4x4 complex mass matrix Y, which reads, 
in the {B, W^, H^, H^) basis: 

/Ml -Mzcpsw 

M2 Mzcpcw 

-Mzcpsw Mzcpcw 

\ MzSjsSw -Mzspcw -A* 

where Mi is the bino mass parameter, Mz is the mass of the Z boson and sw = sin 9w is the 
sine of the weak mixing angle. We adopt the abbreviations = cos (3 and sp = sin /?. Due 




Mz-spsw \ 
-Mzspcw 
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(5) 
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to the Majorana nature of neutralinos, only one 4x4 complex unitary matrix N is required 
to diagonalise Y, where N*YN'' = diag(m^o, m^o, m^^o, m^o). Besides parameters from other 
sectors, the masses and mixings of neutrahnos and charginos can thus be described by three 
independent input parameters, Mi, M2 and /i. If all three parameters are real, then X and 
Y can also be chosen to be real, while each of the rows of N can be chosen to be purely real 
or purely imaginary such that all neutralino masses are positive. 



2.3 Sfermion Sector 

At tree level, the physical squark and charged slepton states, fi. /2, are the mass eigenstates 
of a 2 X 2 complex mass matrix, which reads in the {fi, Jr) basis for each flavour, 

_ / M| + + M| cos 2p{li - Qf-s^,) rrifX} 
f I rufXf M?^ +m| + M|cos2/3g/4 

with 

Xf = Af - fi* {cot P, tan P} , (7) 

where {cot /9, tan /3} applies for up- and down-type sfermions, respectively. The soft SUSY- 
breaking parameters introduced in the sfermion sector are M| and M^^, which are real, and 
the trilinear coupling Af, which can be complex. The phase (pAf can play an important role 
in loops involving the supersymmetric partners of the heavy third-generation SM fermions, 
t, b, T, where the term mfXf appears in couplings of sfermions to Higgs bosons. 



2.4 Tree- level Decay Width 

For the interaction of neutralinos with neutral Higgs bosons, the relevant piece of the La- 
grangian can be written in terms of tree-level mass eigenstates as, 

^ = ^ hl^^i^RCSno + ^l(-1)^-(-1)'^-CJ,o)] x?, (8) 

where cor/l — |(1 ± 75), and k labels neutral Higgs bosons, i.e. — {h, H, A, G}. A minus 
sign appears between the ujr and ujl terms for the CP-odd Higgs states. The couphngs, 
C^(o , are given by 

where 

Cijhl = [{akNis + hkNii) {swNji - cwNji) + {akNjs + bkNj^) {swNn - cwNi2)] 

ttk = {-Sa,Ca,iSf}^,-iC(3j 

bk = {-Ca, -Sa,-iC/3„,-is^„}. (10) 

The tree-level decay width F*''^^ for the two-body decay x° — > x°/i°, where /i° = {h,H,A}, 
can then be written as 

F*- = Y^\C5,of «:(m|o, m|o, ^0 + m|o - m^o + 2(-l)^-m,om,o], (11) 
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with 



(12) 



In order to obtain a prediction for the decay width at one-loop level, the parameters ap- 
pearing in the lowest-order result and the fields of x?; Xj'; ^fc need to be renormahsed. We 
describe their renormalisation in the next section. Note that the mixing matrix elements 
involving a, and are not renormahsed in our scheme, and /3„ is set equal to /5 only 
after the renormalisation has been carried out. 

3 One-loop Calculation for Xj^a cind Combination 

with Higher-order Contributions 

We have calculated the full one-loop vertex corrections to the process Xj^a, where 

= {^1; ^2; ^s}; taking luto accouut all sectors of the MSSM and the full phase dependence 
of the CP- violating parameters and M3. We assume a unit CKM matrix. Examples 
of genuine one-particle irreducible (IPI) vertex diagrams are shown in Fig. [TK,b,c. Fig. [T]i 
shows an example of a reducible diagram, where a Higgs boson mixes with a Z boson or a 
Goldstone boson. For our calculations we have made use of the program FeynArts, allowing 
automated generation of the Feynman diagrams and amplitudes [28-30]. In conjunction, 
we utilised the packages FormCalc and LoopTools for the calculation of matrix elements 
and loop integrals [31]. For regularisation we use dimensional reduction, according to the 
prescription of Refs. [31,32]. We supplemented the model files available in FeynArts with 
counterterms for the 2- and 3-point vertices involved, specified according to the renormali- 
sation prescription outlined below. 




Figure 1: (a) Example of a IPI vertex diagram for X2 ^ Xi''- involving ferniions and sfermions. There are 
UV-finite subsets of diagrams with the same generations and/or flavours; (b,c) Examples of IPI vertex dia- 
grams involving gauge bosons, Higgs bosons and their superpartners; (d) Examples of reducible G-Z mixing 
self-energy diagrams. The particles are labelled according to / = g, i^, x = x°, x^l^, = /i, iJ, A, G, H^,G^ 
and V = -f,Z, W^. 



3.1 Renormalisation in the Chargino and Neutralino Sector 

A significant number of one-loop calculations have been carried out in the chargino-neutralino 
sector of the CP-conserving MSSM with real parameters, see e.g. Refs. [33-40], with the 
renormalisation schemes of Refs. [39,41] also applicable for complex parameters. More re- 
cently, CP-odd observables have been calculated at one-loop level in the CP-violating MSSM, 
see e.g. Refs. [41-43], but these calculations did not always require a dedicated renormalisa- 
tion scheme as the specific observables calculated were UV-finite. In order to renormalise the 
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fields and parameters in this sector we introduce counterterms and renormalisation constants 
of a similar form to Ref. [37]. However, we apply different on-shell conditions for the mass 
parameters and we extend the formalism to the general case including CT'- violation, prop- 
erly taking into account imaginary parts arising both from the complex MSSM parameters 
and from absorptive parts of loop integrals. In this work, we allow the sfermion trilinear 
couplings, Af, and gluino mass parameter, M3, to have CP-violating phases, as inspired by 
the CPX scenario. The full complex MSSM also contains phases in the ncutralino sector at 
tree level, where Mi, M2 and may be complex (there are only two physical phases, since 
one of the phases of Mi and M2 can be rotated away) . For the numerical results obtained in 
this paper we do not need to specify a renormalisation scheme for these phases, since they 
are zero in the CPX scenario. A discussion of this issue will be deferred to a forthcoming 
publication. The mass matrices each receive a counterterm as follows. 



X ^X + 5X, Y ^Y + 5Y, 



(13) 



where 6X and 6Y are 2x2 and 4x4 matrices, respectively. Their elements contain three 
new renormalisation constants, 5 Mi, 5M2 and 5/1, as well as renormalisation constants from 
other sectors. We introduce renormahsation constants separately for the left and right- 
handed components of the incoming and outgoing fermion fields, as follows. 



\6Z^)ijUjRXj , 



Xi 



Xi 
Xi 



^SZt)ijLOR, 
\^Z^)ijUJLi 



X',UR^X'^{l + \5Z^)^jUJR, 



[I + \5Z^)i^URxl 



(14) 



where j is summed over 1,2 (1,2,3,4) for the charginos (neutralinos) . Note that we have in- 
troduced barred renormalisation constants for outgoing fermions and incoming antifermions. 

In the CP-conserving MSSM, these are related to the non-barred renormalisation constants 
for incoming fermions and outgoing antifermions by a Hcrmitian conjugate, i.e. 6Zij = SZjj. 
For the CP-violating MSSM, we choose to treat these quantities as independent at this stage, 
with more discussion to follow. 

Inserting the above transformations into the Born Lagrangian, we obtain formulae for 
the renormalised IPI two-point vertex functions, Fij = iS~j^, where Sij is the loop- corrected 
propagator and hatted quantities arc renormalised. The field renormalisation constants arc 
then fixed by requiring that these i^j are diagonal for on-shell external particle momenta 
and that the propagators have unity residues, namely. 



lim 



rijXi{p)\p2=rnl = 0, 
'■^iiXiip) Xii 



Xi{p)ri 

hm htXiip) 



Xi 



0, 

Xii 



(15) 
(16) 



where Xi = Xi {hJ = 1)2) or Xi {h3 = 1)2,3,4) and i 7^ j. We also require that the 
loop-corrected propagator has the correct Lorentz structure in the on-shell limit. Namely, if 
we decompose the self-energies as 



^^AP') =/^lS{;.(/) +/a;HSg(p2) +a;^Sf (/) (p^) 



(17) 
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then, in order to ensure Su has only a scalar and vector part on-shell, we require, 

Eij(m|J = E^KJ, (18) 
K) = EfK). (19) 



The conditions in Eqs. (jlSl) . (1161) and ( 1191) lead to the following off-diagonal and diagonal 
chargino field renormalisation constants, respectively (Eq. (fTSl) is then automatically satis- 
fied). 



= 2- ["^^±S^''/^(mL) + m~±m-.±J:^{^Jml±) + m-±E^^/"^^(mL) 

+m,.E!^f ^mj.) - m^.^ {U*6XV^)^^ - m-.^ {V5XW)^;\ , (20) 

Xj Xi 

+m^^E'J}f\ml^)-m-^^iU*SXV%-m^^^iV^^ (21) 
SZ^^u = -^f^irriL^-rriL. [E^;,(m^.) + E«;,(m^.)] -m,. [E!:<;(m^.) + E!5(m^.)] 
■ ^ [E!;;,(m|±)-E^^,(m|±) + (r5Xtf/^),,-(f/*5X\/t),,], (22) 



'-ff = -^-^(-^-)— ..K-) + S-:.K-)] [S!5(m^,) + E!5(m^, 

T^[E!f;,(m^±)-E^5,(m^±) + (l^5Xtt/^),,- ([7*5X0)^^] , (23) 

Xi 



while the field renormalisation constants for neutralinos are given by 

+^x?So,S^''(^o) - m,o^^ (X*5rxt)^^. - m,o^^ (iV^yt^r^) J , (24) 
'^^f = ^^S" = -S5:{f (m|o)-m|o [E£,(m|o) + E£,(m|o)] -m,o [EJ^; (m|o) + Ej^'(m|o)] 

±^KiKo)-S^S(m|o) + (X5FtAr^),,-(X*5rxt),,]. (25) 

Xi 

Here E^^(m^) denotes the derivative ^^^^-^ |p2=^2 , and E_ and Eq indicate chargino and 
neutralino self-energies, respectively. The relations between the left- and right-handed con- 
stants for the Majorana neutralinos result from charge conjugation symmetry. For both 
the charginos and neutralinos, the barred constants, bZ^j^ ^ differ from [SZ^J^Y in their 
absorptive parts only. In the CP-conserving MSSM, this difference vanishes. Also, the 
terms in the diagonal constants vanish if there are no complex parameters, and we recover 
the formulae from Ref. [37] in this case. In the CP-violating MSSM, this term is non-zero 
(this term also appears as a purely imaginary contribution in Ref. [44], with which our 
results agree up to absorptive parts). In the CP- violating MSSM, the absorptive parts of 
loop integrals for unstable particles may enter the squared matrix element at the one-loop 
level since they can be multiplied by imaginary coefficients arising from the complex param- 
eters. In the literature, the issue of the treatment of absorptive parts of loop integrals in 
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field renormalisation constants has found considerable attention, mostly in the context of 
the renormalisation of the SM, see e.g. Refs. [45-50]. A possibility that has been advocated 
for instance in Refs. [45,48] is to discard the absorptive parts of loop integrals in the field 
renormalisation constants, while keeping any complex parameters in the coefficients, indi- 
cated by inserting the symbol Re into the renormalisation conditions in Eqs. ( |T5|) . ( fT6i) . ( lT8i) 
and (fT9l) . With this choice the hermiticity relation, 5Z^-^ = {6Z^J^y , is restored, but one 
must include all reducible self-energy diagrams and may have to introduce additional finite 
normalisation constants to ensure the external particles have the correct on-shell proper- 
ties. Renormalisation conditions without Re were used in Ref. [47] for the SM, as a way of 
ensuring the correct on-shell conditions and gauge-independent matrix elements. Although 
the hermiticity relation between renormalisation constants is not valid in this case, the au- 
thors of Ref. [47] showed that the CPT theorem still holds. Nevertheless, the issue of an 
appropriate field renormalisation of unstable particles on external legs remains under debate 
in the literature. For the class of processes considered in this paper, it turns out that all 
absorptive parts of external neutralino self-energy diagrams cancel when the squared matrix 
element is summed over all spins. This is due to the relation between left- and right-handed 
components of the (Majorana) neutralinos, and does not apply for (Dirac) charginos nor 
for spin-dependent calculations. Hence, for the numerical results presented in this paper 
absorptive parts of loop integrals do not contribute. 

It should be noted that the prescription for the field renormalisation constants given 
above is valid for the most general case of CP-violating parameters in the complex MSSM. 
As mentioned above, we restrict the analyses in this paper to cases where the parameters Mi, 
M2 and /i are real. Therefore we do not specify the renormalisation of CP-violating phases 
of those parameters in what follows below. This issue will be addressed in a forthcoming 
pubhcation. 

For the parameter renormalisation of Mi, M2, /i, we use an on-shell approach, because 
this is convenient in processes with external charginos and neutralinos. In the chargino- 
neutralino sector, we have three independent input parameters. Mi, M2, n, which determine 
the tree-level masses, m^^-, of the six fields, xf,2^ Xi,2,3,4- The loop-corrected masses, M^., are 
then defined as the real parts of the poles of the corresponding loop-corrected propagators 
Sii- At one-loop order they may be written in terms of the renormalised self-energies as 
follows, 

M^, = m^^Jl - ^Re[S,^(my + ±j!iml)]] - iRe[Ef (my + S^KJ]- (26) 

We fix three of these masses on-shell by requiring that the pole masses M^^ coincide with 
their tree level values my.. This gives us three equations to solve for 6 Mi, 6M2 and 6fi. 
The remaining three masses will be different to their tree-level values. There is obviously a 
freedom of choice here in the three masses that are used in the on-shell conditions. It should 
be noted that the "most convenient" choice for those masses will depend on the process 
under consideration and may even be different in different regions of parameter space. In 
Ref. [37], the masses of Xi, xf and xf were fixed on-shell. This choice is advantageous for 
processes where charginos appear as external particles, and ensures a proper cancellation 
of the infra-red divergences present in QED corrections. For the processes considered in 
the present paper, it is convenient to have the two lightest neutralinos on-shell. We thus 
choose to fix the masses of Xi, X2 and xt on-shell. We found this to give numerically stable 
results for the hierarchy of Mi < M2 <^ /i among the mass parameters in the chargino- 
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neutralino sector, while for other processes and parameters we found that different choices 
can be favourable. The resulting expressions for 5Mi, 5M2 and 5^ are given below, 

5Mi = [2(iVi3iVl4iV22 - iVl2A^23iV24)q2) + V^22iV22 + 2t/2ll^2liV23iV24)iV(l) 

-(f/22V^22iVi'2 + 2f/2lV2liVl3iVl4)iV(2)]/i^, (27) 

5M2 = [2{Nl^N2^N2A-N^^Nr^Nl^)C(2)~U22V22NlN^^^^-U22V22NlN(2)\/K, (28) 
5/i = [-{Nl2Nl^-NlNl2)Ci^2) + U2iV2iNl^N(^-^-U2iV2iNlN^2)\/K, (29) 

where 

q,) = Re[m-±[S,^(m2±) + Sf(m^±)] + Sf (m2±) + S,f (m^i)] 

-25X2iU,2Va-25X^2UaVi2, 
iV(,) = ReKo[S,^(m|o) + S,^(m|o)] + Sf (m|o) + S,f (m|o)] 

-4<5Fi3iV^iiV.3 - 45F23iV.2iV^3 - 45ri4iV,iiV,4 - 45F24iV,2iV,4, 
= 2[/22V22(A^iiAr| - Nl2Nl^) + 4f/2iV2i(A^iiAr23Ar24 - Ni^Ni^Nj^). (30) 

Note that, since we are assuming Mi, M2 and /i to be real, the mixing matrix elements 
always appear in combinations where the conjugate is not needed. 



3.2 Renormalisation in the Higgs Sector 

For the Higgs sector we follow the renormalisation scheme of Ref. [20]. The independent 
parameters of the Higgs sector are taken to be Mh± and tan/?. One field renormalisation 
constant is introduced for each Higgs doublet, and tan/3 receives a counterterm as follows 

^1,2^(1 + ^5^^1.2)^1,2, tan^ ^tan/5(l + 5tan/5). (31) 

As in Ref. [20] we adopt DR renormalisation for the fields and tan (3, where the counterterm 
for the latter is given by 5tan/5°^ = ^{5Z^ — 5Z^), and on-shell renormalisation of the 
charged Higgs boson mass, 

5Ml± = Re Eh+h- {Ml± ). (32) 

The loop-corrected neutral masses Mh^ are then defined as the real parts of the poles of the 
diagonal elements of the 3x3 Higgs propagator matrix, as in Ref. [20], with M^^ < Mh^ < 

The correct on-shell properties of Higgs bosons appearing as external particles in physical 
processes, and thus a properly normalised S-matrix, are ensured by the introduction of finite 
wavefunction normalisation factors Zij. These Z-factors are a convenient way to account for 
the mixing between the Higgs bosons and to incorporate leading higher-order contributions. 
Following Ref. [20], a renormalised IPI vertex Fi with an external Higgs boson i {i = h, H, A) 
then takes the form 

\J~^i{ri + ZijFj + ZikFk + ...). (33) 

Here j, k are the remaining two of /i, if, A and are not summed over, and the ellipsis refers to 
the mixing contributions with the Goldstone and Z bosons which we will consider in Section 
13.41 Without the ellipsis, the normalisation of the wavefunctions can be expressed in terms 
of a 3 X 3 non-unitary matrix Z, where Zjj = \fZiZij and Zjj = 1. We use the same formulae 
for the Z matrix elements in terms of the Higgs self-energies as those derived in Ref. [5]. 
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3.3 Renormalisation in Other Sectors 



We parameterise the electric charge, e = V4vra, in terms of a{Mz) = a(0)/(l — Aa), where 
Aa = Attiept + Aa^J^ is the shift in the fine-structure constant arising from large logarithms 
of light fermions. This yields the following counterterm for the charge renormalisation, 

- 2n,(0) - 2 • (34) 

Here n^(0) = ^^^p^|fc2=o is the photon vacuum polarisation, and the large logarithms 
involving light fermion masses drop out in Eq. (1341) . The renormalisation constants in the 
gauge boson sector are defined as follows, 

M|^M| + 5M|, ^ M^ + 6M^, Sw^sw + Ssw, (35) 

where 

f6Ml 5Ml 



The renormalisation constants are then deduced from on-shell conditions for the masses of 
the W and Z bosons, 

8M^ = Re S^i^(M^), 5M| = Re Si^(M|). (37) 
3.4 Vertex Renormalisation 

The 3-point vertex for XiXj^fc) where = {h, H, A, G}, can be renormalised by a coupling 
counterterm as follows, 

= ^— +C^^io(5^e-^-^) + -y^(5^f"a^/^ + 5^^^^^ 
2cwsw sw cw 'J'^" 

+\{6Z,o,Cf/,' + SZ.o^Cf/^^ + 6Z,o^Cf/J^ + 6Z,oaCf/Gh (38) 

where 

Scijho = [{akNi3 + hkNii){5swNji - ScwNj2) + {akNjs + bkNj4){6swNii - 6cwNi2)]. (39) 

The 3-point vertex for x^x^ha is then constructed using the 3 x 3 Z matrix for the normali- 
sation of wavefunctions as in Eq. flHHj) . This automatically includes the reducible self-energy 
diagrams involving h, H, A. For a complete one-loop result, reducible diagrams involving 
mixing self-energies of Higgs bosons with the G and Z bosons, such as those in Figure [T](d), 
must also be included. In order to ensure a proper cancellation of the gauge parameter 
dependence, we follow the approach of Ref. [5] and evaluate these reducible contributions, 
pcz.se^ strictly at the one-loop level. Our full result, pf^iiLoop tj^en obtained by combining 
these contributions with those of genuine vertex type, F^^^, as follows, 

where = {h,H,A} are the tree-level states with tree-level masses, m/^o, and are summed 
over. In contrast, M/j^ is the loop-corrected mass of the Higgs boson ha in the physical 
process, i.e. one of /ii,/i2;^3- Numerically, inclusion of the G-Z mixing did not have a 
significant effect. Across the CPX parameter space studied, the effect of this correction on 
the decay widths was less than 0.1%. 
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3.5 Combination with Higher-order Results 

As Higgs propagator-type corrections are known to be large, we have combined our one-loop 
result for the genuine vertex contribution with state-of-the-art two-loop propagator-type 
corrections obtained within the Feynman diagrammatic approach, as implemented in the 
program FeynHiggs [20,22-24]. These contributions incorporate in particular the full phase 
dependence at 0{atas), while we do not include here further two-loop corrections that are 
known only for the case of real MSSM parameters. Using Eq. fl40p . we combine the two-loop 
Z factors and Higgs masses M/j^ from FeynHiggs 2.6.5, with our own genuine vertex (Z"^^^) 
and G-Z mixing (^/"G.z.se-j corrections to the process Xi ~^ X^ha, thereby obtaining the most 
precise predictions for the corresponding decay widths and branching ratios in the MSSM 
with complex parameters. 

In order to investigate the effects of the genuine vertex contributions for the process ~^ 
X^ha we will in the following compare our full result with an Improved Born approximation. 
The latter is obtained by summing over the tree-level amplitudes for x° X^h^y weighted 
by the appropriate Z factors and evaluated at the loop-corrected Higgs masses. 

We will always compare our numerical results to this Improved Born approximation, rather 
than to the strict tree- level result of Eq. ( fTTl) . This allows us to separate out the effect of our 
new genuine (process-specific) vertex corrections from those corrections coming from mixing 
effects and mass shifts in the Higgs sector which are already known to be large. Thus when 
we speak of the percentage effect of our genuine vertex loop calculations on the partial decay 
width, r, we are referring to the ratio 

^ Full Loop Improved Born 

Tlmprovcd Born 

As well as our full MSSM calculation, we will show approximations, where only some (UV- 
finite) sets of diagrams such as third generation quarks and squarks, i.e. t, t, b, b, are included 
in the genuine vertex corrections. In all cases, the two-loop propagator-type corrections 
from FeynHiggs are evaluated in the full MSSM. Various other approximations exist in the 
literature. In Ref. [18], only the one-loop 3rd generation (s) quark contributions in the real 
MSSM were considered. In Ref. [17], all one-loop (s)fermion contributions in the real MSSM 
were considered. Our full results thus go beyond these works, as we include all possible 
MSSM particles in the loops, we allow complex trilinear coupling and gluino parameters, 
and we incorporate complete one-loop and leading two-loop contributions from the Higgs 
sector. 



(41) 



(42) 



4 Numerical Results for the Decay Width 

As explained above, we will discuss in particular the case of the CPX benchmark scenario [4], 
which gives rise to an unexcluded parameter region with a light Higgs. We use the following 
parameters for the CPX scenario unless specified otherwise, 

CPX: /i = 2 TeV, MgusY = 500 GeV, IM3I = 1 TeV, \Af \ = 900 GeV, 

0Af3 = 0A.b..=^/2, M2 = 200 GeV, M, = (5/3)4M2, = 172.4 GeV, (43) 
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X3,4iX2 






Mr 

tl,2 


Mr 

Ol,2 






~2001.2(-5),2003.1(-2),2003.4(1) 


1000 


-500 


332,667 


471,531 


198.5,198.5 


94.7 



Table 1: Masses in GeV of sparticles in the CPX scenario with tan/3 = 5.5, where the neutrahno and 
chargino masses are the tree-level values used in the loops. The number in brackets is the loop-correction to 
the last digit, evaluated from Eq. (pS)) . 



where Msusy = Ml = Mj^, see Eq. (j6]), and tw = swjc-w- The values given in Eq. (H3l) 
differ from the ones given in Ref . [4] in the value of the top-quark mass and in value of | | , for 
which an on-shell value is used that is slightly shifted from the DR value specified in Ref. [4] 
(see also Ref. [5]). In Eq. fHHj) we specify MgusY and \ Aj\ for all sfermions, although it is only 
the third generation which plays a significant role in Higgs phenomenology. For sfermions 
in the first and second generations, EDM constraints are more stringent (see Ref. [51] for a 
recent study and references therein), so we set the corresponding phases of to zero in the 
first two generations. The value of M2 does not play a large role in Higgs phenomenology; 
we choose the nominal value of Mi = 200 GeV in order to agree with other studies, but we 
will investigate how its variation affects our results. We will also show the effect of varying /i 
and At,b,T- The masses of the supersymmetric particles in the CPX scenario of Eq. ( l43l) are 
given in Table [H For comparison, we also present numerical results for the CP-conserving 
case. In particular, we consider the "small aes scenario" [52], with the following parameter 
values 

small OefT : /i = 2 TeV, MgusY = 800 GeV, IM3I = 500 GeV, Xf = -1.1 TeV 

M2 = 500 GeV, Ml = (5/3)t^M2, = 172.4 GeV. (44) 

The small aes scenario has some similarities to the CPX scenario, including large values of 
fi and \Af\, where Aj is related to Xf according to Eq. (JTj). As usual for these benchmark 
scenarios, the parameters characterising the Higgs sector at lowest order, i.e. Mh± (M^) and 
tan/? in the case of the CPX (small Oeff) scenario, are varied. 

We furthermore investigate a specific case of a CP-conserving scenario giving rise to a 
very light Xi^ inspired by a recent study [53] which showed that very light neutralinos are 
not ruled out by experimental data. Here the GUT relation between Mi and M2 is relaxed, 
allowing Mi to be chosen to be such that the lightest neutralino is approximately massless. 
Unless stated otherwise we use the following parameters, 

light Xi- = 600 GeV, Mgusv = 500 GeV, IM3I = 1 TeV, tan/3 = 20 

Af = 1 TeV, M2 = 400 GeV, m^o = 0, = 172.4 GeV. (45) 

We start with numerical results for our genuine vertex corrections to the decay width for 
X2 — >■ Xi^i in the CPX scenario. Figure 2(a) shows the partial decay width r(x2 ~^ Xi^i) 



as a function of M/^^. The value of tan/5 is fixed at 5.5, while Mjj± is varied as input. 
The dotted Improved Born curve shows the result obtained by combining the tree-level 
amplitudes with 2-loop Z matrix elements and masses according to Eq. (HTl) . The other 
curves incorporate our new results for the genuine vertex corrections, taking into account 



different sets of loop contributions. Figure 2(b) shows the ratio r, defined in Eq. f l42l) . of the 



genuine vertex corrections relative to the Improved Born result as a function of M^-^. We 
see from the figure that the impact of the genuine vertex corrections on the decay width 
is very large. The corrections from the full MSSM contributions to the vertex amount to 
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about 45% for Higgs mass values in the region of the "CPX hole", i.e. for Mh^ ~ 40 GeV. As 
expected, the dominant effect arises from the triangle diagrams containing third generation 
quarks and squarks {t,i,b,b), due to the large top Yukawa coupling, yielding a correction 
of about 35% compared to the Improved Born result. The other (s)fermions also play a 
non-negligible role, in particular through their couplings to neutralinos, increasing the total 
(s)fermion contribution to just under 50%. The vertex corrections from the remainder of the 
particles in the MSSM, namely the vector bosons, Higgs bosons, neutralinos and charginos, 
are negative and contribute about a 5% correction. A similar pattern of the relative impact 
of the corrections is observed if tan f3 is varied while Mh± is adjusted to keep M^^ constant 
at 40 GeV (not shown in the plot). Values of tan/5 below 5 yield a significant increase of the 
decay width. 
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Figure 2: Results for (a) the decay width r(x2 ^ Xi^i) and (b) the ratio r = (Fpuii Loop - 
rimprovcd Born ) / rimprovcd Born in the CPX scenario plotted against Mh-^ for tan/3 = 5.5. {Mu± was var- 
ied as input.) The different curves indicate the inclusion of various subsets of diagrams. 



Such large effects from the genuine vertex corrections are not unexpected in the CPX 
scenario (see also Ref. [5] for an analysis of genuine vertex corrections to Higgs cascade 
decays). It is well known that loop corrections in the Higgs sector can be large, especially in 
this rather extreme scenario with large trilinear couplings and CP-violating phases. Such a 
large value of /x also enhances the effect of loop corrections in the neutralino sector. In Figure 



3(a) we see how the effect of the genuine vertex corrections is further enhanced to values of 
60% or more if /i is increased compared to its value in the CPX scenario of /i = 2 TeV. On 
the other hand, if /i is decreased one obtains correspondingly smaller corrections. 

We also examined the effect of varying the magnitude and CP-violating phase of the 
trilinear coupling. At - 
plot r as a function of 
\At\ = 900 GeV. At 



Ab = A-r, for the third generation of sfermions. In Figure 3(b) 
^At for various values of \At 



we 



First we discuss the bold curve, where 
^At = 7r/2, the loop corrections show a steep dependence on the phase, 
(pAt, emphasising the importance of including these phases in the calculation. At this value, 
hi has its largest CP-odd content, i.e. |Zi3| is largest, while the CP-even contributions are 
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Figure 3: The ratio r = (rpun Loop - rimprovod Born )/rimprovcd Born for X2 ^ Xi^i in the CPX scenario. 
(a) r plotted against /i, with various subsets of diagrams included. Mh± was adjusted in order to keep 
= 40GeV constant and tan/3 = 5.5. (b) r plotted against 4>Af The solid, dashed and dotted curves 
correspond to \At\ = 900, 700, 500 GeV, respectively. Mh± was adjusted in order to keep M^^ = 45 GeV 
constant and tan/3 = 7. (A Higgs mass of Mh^ = 40 GeV was not theoretically accessible for all (pAt when 
\At \ = 500 GeV.) 



suppressed, giving rise to corrections of order r ~ 45% (see Figure 2(b)). When = 0, 
the loop corrections are found to be somewhat smaller, leading to r ~ 30%. On the other 
hand, the effect of the genuine vertex corrections is maximised for (p^^ = vr, i.e. At = —\At\. 
This corresponds to a maximum in |Zii|, so that the lightest Higgs boson is mostly CP-even. 
The genuine vertex corrections for a CP-even Higgs are larger than for a CP-odd one, so 
that their effect is maximised here. Hence the corrections in such a CP-conserving scenario 
can even exceed the ones in the CPX scenario. It should be noted in this context, however, 
that such a light CP-even Higgs boson is of course experimentally excluded. For smaller 



values of \At\ (dotted curves in Figure 3(b)) the corrections are in general smaller, and the 
variation with the phase of At is less pronounced. Nevertheless, even for |At| = 500 GeV we 
find r ~ 35% and r ~ 40% at (pAt = vr/2 and (pAt = tt, respectively. 

We next consider the small a^s scenario. Like the CPX scenario, this scenario has large 
fi and large, negative At. For the small Oefr scenario with Mfj± = 220 GeV and tan/5 = 10, 
we find genuine vertex corrections of size r ~ 35%. The variation with /i, shown in Figure 
results in a pattern that is very similar to the one observed for the CPX scenario in 
The size of the correction scales approximately linearly with /x, and the inclusion of 
the full (s)fermion contributions yields a shift of about 10% compared to the contribution of 
only the third generation (s)quarks. The non-(s)fermionic corrections to the genuine vertex 



4(a; 



3(a, 



give rise to a downward shift of about 5%. In Figure 4(b) the small scenario is modified 
by varying the phase (pAt while keeping \At\ = \Xt — fi* cot f3\ constant. We find that, like 
the CPX scenario, the genuine vertex corrections have the largest effect of order 35% at 
the nominal value of (pAt = tt, while the corrections are only a few percent when the phase 
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is maximally CP- violating for (^At = 7r/2. This can again be compared with Figure 3(b) 
Unlike the CPX scenario for which the corrections are minimised at (pAt = 0? 27r, in the 
small OefT scenario the vertex corrections exhibit another extremum here, with r ~ —20%. 
As for Figure E^b), the dotted curves in Figure 4(b) show the reduced effect of the loop 
corrections when \At\ is decreased. Here we vary tan/? in order to produce the desired \At\ 
from Xt = -llOOGeV using Eq. (El). 




(a) (b) 



Figure 4: The ratio r = (FfuII Loop - Timprovod Born) /rimprovod Born for ^2 Xi^^i in t^c Small ckoff scenarlo. 
(a) r plotted against fi with various subsets of diagrams included. In this plot Mh± = 220 GeV and tan/3 = 
10. b) r plotted against for three different values tan/3 = 3.3,5, 10 and hence \At\ = 500, 700, 900 GeV 
respectively (see Eq. dS}), with Mh± = 220 GeV. 

In addition to the CPX and small aes scenarios, we examined all of the kinematically 
open decay modes of the form Xi X^h for the SPS benchmark points [54]. Corrections of 
over 10% to the partial decay widths were found to be common, indicating that significant 
effects are not limited to scenarios with very large values of fi. Another scenario which we 
examined was the benchmark point LM5, studied in the CMS Technical Design Report, in 
the context of the decay X2 ~^ Xi^ [10]. For this scenario we found the corrections to the 
partial decay width to be around 5%. However, due to the large branching ratio of 85% for 
the process, these corrections translated into an effect of less than a percent on the branching 
ratio. 



5 Phenomenology at the LHC 

5.1 Numerical Results for Branching Ratio 

In the previous section, we found that the genuine vertex corrections to the partial decay 
width r(x2 Xi^i) were of order 45% in the CPX scenario. For phenomenology at the 
LHC it is important to consider, in addition to the decay widths, also the branching ratios 
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of neutralinos. In this section, we compute the branching ratios of incorporating our 
loop-corrected decay widths for X2 ~^ Xi^i,2,3- As well as calculating the genuine vertex cor- 
rections to x° X^ha, we have also calculated the genuine vertex corrections to X^^y 
using a similar procedure to that detailed in the previous sections for the Higgs vertex. We 
incorporate these into the branching ratio calculation too. In the CPX scenario, depending 
on its mass, Xo, can decay via the following decay modes: 

Xl ^ X>u X>2, X>3, xlZ, xlff. /i,2/, /i,2/. (46) 

Where kinematically possible, we calculate the decays X2 ~^ Xi^a, which produce on-shell 
neutral Higgs bosons, as two-body decays, including the genuine vertex corrections as de- 
tailed in the previous sections. Where kinematically possible, we also calculate the decay 
X2 ~* Xi^ into an on-shell Z boson as a two-body decay, including the equivalent genuine 
vertex corrections as mentioned above. Finally, we calculate the 3-body decay X2 Xiff- 
For this, we include, firstly, the diagrams where an off-shell Higgs boson is exchanged (i.e. 
where some or all of hi, h2, are too heavy to be produced on-shell). For these diagrams we 
use the unitary U matrix elements and masses from FeynHiggs to construct effective cou- 
plings (see Ref. [20]) which take into account the two- loop Higgs propagator-type corrections. 
Secondly, in the three-body decay, where the kinematics do not permit an on-shell Z boson, 
we include the diagram where a Z boson is exchanged, along with the diagram where the 
would-be Goldstone boson, G, is exchanged (in this way a proper cancellation of the gauge 
dependence is ensured). Thirdly, we include in the three-body decay the diagrams where a 
sfermion is exchanged. As the neutralino mass approaches the scale of the sfermion masses, 
the possibility of on-shell production of sfermions arises, which subsequently decay into Xi- 
To take account of this threshold region, we include a finite width for each sfermion, cal- 
culated from its self-energy. All self-energies and two- and three-body partial decay widths 
were calculated with the help of FeynArts and FormCalc. 

Throughout most of the CPX parameter space, X2 ~* x\Z can proceed as a two-body 
decay. For the CPX scenario, we find that the genuine vertex corrections to this decay width 
can be of order 30%. However, the amplitude is suppressed by several orders of magnitude 
in the CPX scenario, since the Z boson only couples to the higgsino component of each of the 
neutralinos, while the large value of ^ renders x? and X2 mostly bino and wino, respectively. 

The resulting branching ratios of X2 in the CPX scenario are plotted as a function of the 
neutralino mass, fn^o^, in Figure [5|i, with tan (3 = 5.5 and Mh-^ = 40 GeV. Both the Improved 
Born and full MSSM vertex-corrected results are shown. We see that for m^o ^ 190 GeV, 
BR(x2 Xi^i) ~ 100%, and therefore the loop corrections to the X2 ~^ Xi^i partial 
width have negligible effect. As one increases m^o from 190 to 470 GeV, the on-shell decays 
X2 ~^ Xi^2 and X2 ~^ Xi^3 become kinematically allowed. This causes BR(x2 Xi^i) 
to vary from 100% to around 25%. In this region, the three competing decay modes into 
Higgs bosons all receive large vertex corrections of order 50%. However, since these vertex 
corrections have similar structure, their effects tend to cancel each other out, producing an 
effect of only a few percent on the branching ratios. Thus, the Improved Born approximation 
works well in this region. 

The effect of vertex corrections on the branching ratio will be more significant in regions 
of parameter space where there is another competing decay mode of X2 which does not have 
loop corrections of a similar structure to X2 ~^ Xi^i- In the CPX scenario, this competition 
will never be provided by the highly suppressed decay into a Z boson. However, for m^o 
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Figure 5: Branching ratio for each of X2 Xi^i,2,3 and for the other decay modes, X2 ~^ Xi^ and 
X2 ~* Xiff (labeUed "Other"); (a) shown as a function of m^o, for Mh^ = 40GeV and tan/3 = 5.5 (M2 was 
varied as input to produce the change in m^o); and, (b) shown as a function of for tan/3 = 5.5 and 

M2 = 200 GeV {M^^ was varied as input). In both plots we show the Improved Born approximation as a 
dashed line and the full MSSM result as a solid line. 



large enough, decays via sfermions become important. While the Higgs bosons require both 
a non-zero gaugino and higgsino component to couple to neutralinos, sfermions couple only 
to the gaugino part. Thus, if sfermion decays can proceed on-shell, they will, in this scenario, 
dominate over the Higgs decay modes, rendering BR(x2 Xi^i) ~ 0- Before this, there 
will be a threshold region, which can be seen in Figure 5(a) for 450 ^ m^o ^ 520 GeV. 
Within this region, the existence of competing decay modes means that the genuine vertex 
corrections are very important. The maximum effect occurs near m^o ~ 500GeV, where the 
positive vertex corrections to the Higgs decay widths result in a reduction of the branching 
ratio BR(x2 Xiff) of more than 10% compared to its Improved Born value. 



In Figure 5(b) , we show the branching ratios of X2 a function of M/i^, to be compared 
with Figure 2(a) Here M2 = 200 GeV and tan/? = 5.5, so only decays into Xi^i; Xi^2, Xi^ 
and Xiff (the latter two labelled "Other") are kinematically open. For Mh^ > 50GeV, the 
second lightest Higgs boson is too heavy to be produced on-shell and so BR(x2 Xi^i) is 
close to 100%. In the CPX hole, with Mh, ^ 40GeV, we find BR{x2 ^ X^h) ~ 79%, an 
increase of around 3% compared to the Improved Born value. 

Thus, although we found large loop corrections to the partial decay widths of X2 ~* 
Xi^i,2,3 in the CPX scenario, the effects on the branching ratios turn out to be significantly 
smaller, because the decays are not competing with other modes and so the large genuine 
vertex corrections cancel each other out. This will also be the case for the small scenario, 
in which the Z decay mode is also suppressed and the sfermions are heavy. However, this 
situation is not generic, and large vertex corrections can affect the branching ratios if there are 
other competing decay modes with vertex corrections of a different structure. In non-gaugino- 
like scenarios, without a large hierarchy between M2 and fi, the decays into Higgs bosons are 
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more likely to compete with the decays into Z bosons and sfermions. For example, in Figure 
[6] we show the Improved Born and full MSSM branching ratios for the "light Xi scenario" 
of Eqn. US Here we can have BR(x2 Xi^) ~ BR(x2 Xi^)- We computed genuine 
vertex corrections to both T{x2 Xi^) and T{x2 —>■ x\Z)i ^^id found the former (latter) to 
be negative and of order 20% (3%) and 35% (2%) for At = 500 GeV and At = 1200 GeV, 
respectively. The corrections are further enhanced at large values of tan (3. In this scenario, 
the branching ratio for X2 ~^ Xih happens to be near 50%. Thus, the effect of the vertex 
corrections on the branching ratios is maximised in this case. The plots in Figure El show 
corrections to the branching ratio of more than 10%. 
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Figure 6: Branching ratios for each of X2 — > Xi^j X2 ~^ Xi^ and X2 ~^ Xiff (labelled "Other") in the 
CP-conserving "light scenario": (a) shown as a function of tan /3, for fixed Af = ITeV; and, (b) shown 
as a function of At, for fixed tan/3 = 20. We show the Improved Born approximation (I.B.) as the dotted 
line, and the fuU MSSM result (F.M.) as the solid line. 



5.2 Prospects for the "CPX Hole" 

In the previous section, we found that X2 ~^ Xi^i has a large branching ratio, BR(x2 — 
Xi^i) ~ 79%, for the "CPX hole", i.e. in the region where a light Higgs is unexcluded by 
present data. We now investigate whether Higgs production in neutralino decays at the LHC 
could help to cover this parameter region. Consider the SUSY cascade decay chain starting 
with a gluino; 

g^ff^ Xlff - X'JIK - x'jfbbir^T-). (47) 

Coloured sparticles like gluinos are expected to be produced in large numbers at the LHC 
provided they are light enough, (see eg. ReL [55] for detailed analyses of SUSY cascade 
decays). These gluinos will decay into lighter coloured sparticles, namely the squarks with 
masses around 500 GeV (see Tabled]). For most squarks, the only way to conserve R-parity 
will be to decay into X2? xf ^^(^ Xi- As shown in the previous section, 79% of the produced 
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X2 will decay into hi in this scenario for a Higgs mass of 40 GeV. The light Higgs boson 
then decays mostly into bb (91%), and also t~^t~ . 

Branching ratios for all parts of the decay chain, except the decays involving Higgs 
bosons, were computed at tree level using FeynArts and FormCalc. For decays involving 
Higgs bosons, such as t2 tiha, we use an Improved Born approximation, as in Eq. HH We 
computed the branching ratio for g qi^2<i for each of g = m, d, c, s, t, b. We found that 61, 
Ml, (^2, C2, Si, S2 all have substantial branching ratios to decay into X2- Summing over the 
various decay modes we found that 17% of all gluinos produced decay via a squark into Xt 

Combining BR(x2 ~^ Xi^i) ~ 79% with BR(g — > x'il'l) ~ 17%, we estimate that around 
13% of the gluinos produced in this scenario will decay into hi. Thus, SUSY cascade decays 
where a light Higgs is produced in the decay of the second-lightest neutralino appear to be 
a promising possibility to cover this problematic parameter region where standard search 
channels may only have small sensitivities. Detailed experimental analyses would be needed 
to determine whether it is indeed possible in such a case to extract a Higgs signal from the 
SM and SUSY backgrounds^ 

6 Conclusions 

We have obtained complete one-loop results for the class of processes involving the decay of 
a neutralino into a neutral Higgs boson plus a lighter neutralino, x^ — ^ x'jha- The genuine 
vertex contributions to the neutralino decay amplitudes have been combined with state- 
of-the-art two-loop propagator-type corrections for the outgoing Higgs boson. Our results 
take into account all sectors of the MSSM and include the full phase dependence of the 
CP- violating parameters Af and M3. 

For the renormalisation in the chargino-neutralino sector, we have worked out an on-shell 
scheme which properly takes into account imaginary parts arising from complex parameters 
and from absorptive parts of loop integrals. In this scheme in- and outgoing fermions re- 
ceive different field renormalisation constants. Since we have concentrated in this paper in 
particular on the CPX benchmark scenario, where only the parameter M3 and the trilinear 
couplings of the third generation fermions are complex, we have not specified the renor- 
malisation of the phases of the parameters appearing in the neutralino and chargino mass 
matrices. This issue will be addressed in a forthcoming publication. 

In the CPX scenario we find corrections to the partial decay width for X2 ~^ Xi^i of about 
45% relative to the Improved Born approximation. These corrections, which characterise the 
impact of the genuine vertex contributions, scale almost linearly with the higgsino parameter, 
fi. We also found a strong dependence of the size of the genuine vertex corrections on 
the absolute value, l^^l, and the CP- violating phase, (pAt, of the third generation sfermion 
trilinear coupling. The corrections turn out to be even larger for the (CP-conserving) case 
where (p^t = ^ compared to the case of the maximally CP- violating phase of (p^t = vr/2. 

^It should be noted in this context that the CMS collaboration has performed a full detector simulation 
and event reconstruction for the production of a Higgs boson at the end of a cascade of supersymmetric 
particles starting with squarks and gluinos [10]. These results, obtained for the benchmark point LM5, 
cannot be directly translated to the case of the CPX scenario, since in the case of LM5 the Higgs boson is 
much heavier, nih ^ 115 GeV, than in the region of the CPX scenario that we are considering here. The b 
jets resulting from the Higgs decay in the CPX scenario are therefore softer than for LM5, so that cuts on 
the energy of the jets will be less efficient to suppress the QCD background. 
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We have also investigated CP-conserving scenarios, in particular the small aes scenario. 
Similar to the CPX scenario, the small scenario is characterised by large values of fi 
and \At\, and we found large corrections to the partial decay width of about 35%. In both 
the CPX and small scenarios, the predominantly gaugino-like character of the two light 
neutralinos reduces the effect of the large vertex corrections on the branching ratios down 
to a few percent. However, we also applied our results to the "light Xi" scenario and found 
corrections to the branching ratio of more than 10%. We also investigated a number of 
other CP- conserving scenarios and found non-negligible corrections. Our results can also be 
applied to the process ha x'iX'j- 

Based on our precise predictions for T{x2 Xi^i) and the corresponding branching ratio, 
we have investigated the prospects of SUSY cascade decays, where a Higgs is produced in the 
decay of a neutralino, for covering the parameter region of the CPX scenario in which a light 
Higgs boson is unexcluded. We find that around 13% of all gluinos produced at the LHC in 
the CPX scenario will decay via the second lightest neutralino into the lightest Higgs boson. 
Thus, Higgs production in neutralino decays looks promising as a search channel for such a 
light Higgs, while standard search channels may have small sensitivities in this parameter 
region. The results obtained in this paper will be provided as a public tool with the aim of 
facilitating further experimental studies of this potentially interesting channel. 
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